#include <iostream>
#include "LSIcompute.h"

int main(int argc, char ** argv)
{
	int m = 10;
	int n = 6;

	double M[] = {
		24, 21, 9, 0, 0, 3,
		32, 10, 5, 0, 3, 0,
		12, 16, 5, 0, 0, 0,
		6, 7, 2, 0, 0, 0,
		43, 31, 20, 0, 3, 0, 
		2, 0, 0, 18, 7, 16, 
		0, 0, 1, 32, 12, 0, 
		3, 0, 0, 22, 4, 2, 
		1, 0, 0, 34, 27, 25,
		6, 0, 0, 17, 4, 23
	};
	
	double * testM;
	int * Groups;

	//testM = new double * [m];

	Groups = new int [m];

	/*for(int i = 0; i < m; i++){
		testM[i] = new double [n];
		for(int j = 0; j < n; j++){
			testM[i][j] = M[i][j];
		}
	}*/
	int pseudoTermNum = 2;
	double dError = 0.08;
	testM = M;

	LSIcompute oLSI(testM, m, n, &dError, Groups, &pseudoTermNum);

	oLSI.runCompution();

	return 0;
}